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An optimized Rayleigh-Schrodinger expansion scheme of solving the functional Schrodinger equa- 
tion with an external source is proposed to calculate the effective potential beyond the Gaussian 
approximation. For a scalar field theory whose potential function has a Fourier representation in 
a sense of tempered distributions, we obtain the effective potential up to the second order, and 
show that the first-order result is just the Gaussian effective potential. Its application to the A</> 4 
field theory yields the same post-Gaussian effective potential as obtained in the functional integral 
formalism. 
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I. INTRODUCTION 



> 

' Since the mid 1980s, the effective potential (EP) in the quantum field theory (QFT) beyond the Gaussian approx- 
.imation (GA) has received much attention because it collects merits and removes weakpoints of the conventional 
'perturbation theory and the GA as shown by the so-called variational perturbation theory in other fields Q|. Sonic 
schemes have proposed to calculate such an EP, for example, designing some non-Gaussian trial wavefunctionals f|](i), 
-using the Brueckner-Goldstone formula with a variational basis (^](ii), rearranging loop diagrams in the functional 
integral formalism with the background field method |Q](iii), optimized expansions in the functional integral formalism 
with the steepest-descent method [Q(iv) and with the background field method Q (v) , and so on. Noting that the 
Q_! 'well- developed Rayleigh-Schrodinger (RS) expansion in quantum mechanics [[| was generalized to QFT Q, recently 
<X) we have developed a variational perturbation scheme with the RS expansion to calculate the EP as per the variational 
^3 'minimum definition jfjj. In the scheme, a free-field Hamiltonian with a mass fixed from the GA was adopted as a 
solvable part so as to perform the RS expansion pi. Obviously, this scheme amounts to a series expansion around 
■the Gaussian EP, and so do those schemes in Ref. |y] (i) — (iii). In the present paper, taking a free-field Hamiltonian 
instead with an arbitrary mass parameter [i as the solvable part, we will apply the RS expansion to solve the functional 
^ Schrodinger equation in presence of an external source so that the vacuum energy functional of the source can be 
obtained. Then, adopting the alternative, but equivalent definition of the EP in the functional Schrodinger picture ||, 
we extract the EP through a Legendre transform of the energy functional. In the new scheme, /i will be determined 
according to the principle of minimal sensitivity JtJ and, consequently, the value of /i in an approximation up to one 
order is different from those up to other orders of the expansion. This way of fixing fi makes the expansion give the 
Gaussian EP with its lowest orders, but be not an expansion around the Gaussian EP. Furthermore, unlike Ref. ||, 
the vacuum expectation value <p of the field operator will naturally be given in the present scheme. Hereafter, we will 
call this scheme the optimized RS expansion (ORSE). Since, for the case of scalar field theory, the EP beyond the 
GA was given only for the A0 4 model except for the <f> 6 model Q(1985) up to now, we use the ORSE in this paper 
to give the EP for a generic class of scalar field models (see Sect. Ill) up to the second order. The resultant formula 
can easily be used to give the EP beyond the GA for a number of concrete scalar field models, for instance, models 
with polynomial or/and exponential interactions. We also show that its application to A0 4 field theory gives the same 
post-Gaussian EP as Ref. g](iv,v), and yields the result in Ref. [|| [||(ii,iii) if one chooses to use the same constraint 
on ip as in Ref. 
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Next, for the sake of convenience, we will first give the complete basis set for a free field theory with an external 
source, which will be employed in the ORSE. In Sect. Ill, the ORSE will be proposed. In Sect. IV, we will perform 
the ORSE to obtain the EP for a class of scalar field theories up to the second order, and its application to the 
A0 4 field theory will be given in Sect. V. Conclusions will be made in Sect. VI with discussions on some possible 
extensions and developments of the present work. 



II. A FREE FIELD THEORY WITH AN EXTERNAL SOURCE 



In this section, we discuss the free-field theory with an external source J x = J{x) in a time-fixed functional 
Schrodinger picture. The Hamiltonian is given by 

= J^nl + \{d x <p x ) 2 + \v?<& - JA X - \f xx + \J^ J x h£j y ], (l) 

where x = (x 1 , x 2 , ■ ■ ■ , x D ) represents a position in D-dimensional space, J = J d D x, fi an arbitrary mass parameter 
and 4> x = 4>{x) the field at x. Tl x = — ijj- is canonically conjugate to <p x with the commutation relation, [0 x ,IIj,] = 
iS(x - y). In Eq.(l), f xy = ( y/-d% + fi 2 )S(x - y) with J z f xz f~J = 5(x - y), and h xy = (-d 2 + fi 2 )6(x - y) with 
j z h xz Ky =S(x-y). 

The functional Schrodinger equation for Eq.(l), i/^^n; J >= E n ° [J]\n; J >, is easily solved || (Here, the 

subscript n in En\j] is the index of eigenstates, and the superscript "(fc)" means "at the kth order of 5". See 
the next section.). First, the eigenenergy of the vacuum state vanishes, and the corresponding wavefunctional is a 
Gaussian-type functional 



|0; J >= 7Vcxp{-i [ (<p x - [ h x ^J z )f xy (cf> y - I h y }j z )} , 

J x.y J z J z 



(2) 



(0) (0) (0) (0) 

where M is the normalization constant (i.e., < J;0|0; J > = 1). One can show that < J;0\(j) x \0; J >= J z h xz J z 
which, unlike Eq.(3) in Ref. Q, is dependent on x. Then, for the above vacuum, the annihilation and creation 
operators can respectively be constructed as 

AfiP; J) = ( 2(27r)W(p) )1/2 L e ' lPX[f{p)[4)x ~ J K * Jz) + lIix] (3) 

and 

At f {P '' J) = ( 2(27r)^/(p) )1/2 / e *™^» ~ / J ») - ^ 

with [Af (p;J),A^(p';J)] = S(p' - p) and A f (p)\0 >= 0. It is not difficult to verify that = 

J d D pf(p)A^(p; J)Af(p; J), where f(p) = ^/p 2 + p? arises from f xy — J d D pf(p)e zp( - x ~ v ^ with p = (p 1 ,p 2 , ■ ■ ■ ,p D ). 
Consequently, the cigcnwavcfunctionals for excited states can be easily written as 

(o) 1 -A- + (0) 
\n;J>=—l[A}(p l ;J)\0;J>, n=l,2,---,oo (5) 

and the corresponding eigenenergies are 

n 

E^[j]=j2f(Piy (6) 

(0) (0) , (0) (0) 

Evidently, the eigenwavefunctionals \n; J > and |0; J > are orthogonal and normalized, < J;m\n;J > = 
{ ran ^i Spifn) Ilfe=i 3(p'k ~ Pin)- Here, Pi(n) represents a permutation of the set {ik} = {1,2, •••,«} and the 

(o) 

summation is over all Pi(n)s. \n;J > describes a n-particle state with the continuous momenta Px-iPi, • • • ,p n - 

(0) (0) , „ (0)(0) 

|0; J > and \n; J > with n — 1,2, • • • ,oo constitute the complete set for Hq, and satisfy the closure |0; J > < 
J; 0| + EZi I d D Pid D p 2 ■ ■ ■ d D p n \n: J >< J; n\ = 1. 



2 



III. OPTIMIZED RAYLEIGH-SCHRODINGER EXPANSION 
FOR THE EFFECTIVE POTENTIAL 



The EP for a field system is equivalently defined through the Feynman graphs, the operator representation, the 
path integral p| , or the minimum expected energy in a set of normalized states j|,|8| . They were used to give the 
loop or Gaussian EP and propose those schemes in Refs. Yet another equivalent definition of the EP is given 

through the vacuum energy functional of an external source obtained by solving the relevant functional Schrodinger 
equation ||. Based on it, we will construct the ORSE in this section. 

In this and next sections, we work with a scalar field model whose Lagrangian density is 



C = -d^ x d^ x - V((f> x ) . (7) 



In Eq.(7), the model potential is assumed to be written as V{4> x ) = J ■^=V{Vl)e in ^ :l , at least, in a sense of tempered 



/2tt 

distributions jL0|. It represents several scalar-field models, such as A0 4 model [p|JI]J|, general and special </> 6 models 
sine-Gordon and sinh-Gordon models |l3), massive and double sine-Gordon model Jl4j], Liouville model JT3j , as 
well as two generic models investigated in Ref. Jl6|] . 

For the system, Eq.(7), the time- independent functional Schrodinger equation in the presence of an external source 
J x is 



(H- { JM\* n >=E n [J}\V n > 

J X 



(8) 



with the Hamiltonian H = / x [^n^ + \{d x <t> x ) 2 + V(<j> x )]. Here, the eigenvalue E n [J] is a functional of J, x . 

For our purpose, Eq.(8) will be modified. We make a shift <j> x ~ * <f> x +$ (3> is a constant), and Eq.(8) can equivalently 
be rewritten as 

[H(4> x + $) - / J x (cf> x + + $, J] = E n [J; $]# n [<fe + J] (9) 

J x 

with H((f) x + $) = ^[ill^ + \{d x cj> x ) 2 + V(cj) x + $)]. This shift is really in the spirit of the background-field method 
|l7j Q(iii). Further, normal-ordering the Hamiltonian in Eq.(9) with respect to a normal-ordering mass M |18|, and 
inserting a vanishing term f x [^H 2 (f> x — ^H 2( f> x ] with /i an arbitrary mass parameter into the Hamiltonian |fl9|] , one can 
have 

N M [H{<j> x + $) - / J x {<j> x + *)] = + H?< * - C (10) 

J X 

with 

H?* = J{-\^ 2 <fl+MM[V{^ x + *)]} (11) 

and 

f 1 1 ( 1 M 2 

C = J [~-f xx + -J J x h~ij y + -I Q (M 2 ) - —h(M 2 ) + J X M . (12) 



Here, the notation Afja[' • •] represents normal-ordered form with respect to M, I n (Q 2 ) = J j + Q?) n , an d 

Mm[V{4> x + $)] = / ^V(n)e in ^+^ +M r I ^ M ^ g. From now on, we will use the normal-ordered Hamiltonian in 
Eq.(10) instead of the original one, which will naturally make the EP in (1+1) dimensions free of explicit ultraviolet 
divergences |9||lS| . Noting that C is a constant independent of <p x and is an exactly-solved Hamiltonian, we can 
formally treat Hj^ as a "perturbed" interaction in the RS expansion To mark the order of the RS expansion, 
an index factor 5 will be inserted in front of Hj ' in Eq.(10). Consequently, Eq.(9) is modified as 

[ftf" + 5H?<*]V n [<f> x + J; 8} = (E n [J; 6} + C)V n [<j> x + $, J; 6} . (13) 

Now, applying the RS expansion, one can solve Eq.(13) to get energy eigenvalues, E n [J;<i>,5], and eigenwavefunc- 
tionals, ^ n [4> x + $>,J;$\. Here, we are interested only in the eigenenergy functional for vacuum state, Eo[J; <i>, $]. 
Obviously, the zeroth-order approximation to Eq[J; $, 5], En [J; <f>], satisfies 
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4°V;0] + c = 4 o) M=o (14) 

and at the nth-order of 8, the correction to Eq [J; <fi] is 

E^[J;$\ = ( <> 0\H^[Q o 1 - (^[J^j-gf*)]"- 1 ^ > (I-".) 



with Q n - T^nId D Pid D p 2 ■ ■ ■ d D Pj \j >< j\. Thus, £ [J; S, 8} = E^[J; 0] + £~ x S"£;< n) [J; 
With the vacuum energy functional So [J; <&, 5], one can have 



= - J v^l [4> x + ®, J-, 8} {4> x + m° 1^ + $, J; S] 



-fx (16) 



where we used the Feynman-Hellmann theorem [Q. Evidently, for 8 = 1, f x =< ^ol^sl^o >• Then, a Legendre 
transformation of Eq[J; <I>, 8] yields the static effective action || 



r e [(p;$,S\=-Eo[J;*,S\- [ Jxfx ■ 

J X 



(17) 



To calculate EP, one can conveniently take f x = 4> in Eq.(16) to fix the arbitrary shifted parameter <&. In analogy 
to Appendix A of Ref. Q (1990) but with the Feynman-Hellmann theorem ^0|, one can show that other choices of 
$ will give rise to the same EP (when the wavefunction renormalization procedure is not needed). Finally, one can 
have the EP f§ 



L 



(18) 

ipx— constant— <1>, 5—1 



If the EP is truncated at some order of 8, the extrapolation of 8 = 1 should be made after renormalizing the 
approximated EP. 

When truncated at a given order of 5, the right side of Eq.(18) will depend upon /i. To obtain an approximated 
EP up to the same order, one can determine /i according to the principle of minimal sensitivity 0. That is, fi should 
be such a value that the approximated EP up to the given order is optimized to be as insensitive to variations in /i 
as possible. This can be realized by analyzing the vanishing first- or higher-order derivatives of the truncated result 
with respect to /i fl S(v). From the next section, one will see that the EP up to the first order is just the Gaussian 
EP. Thus, the ORSE is a systematic tool of improving the Gaussian EP. 



IV. OPTIMIZED EFFECTIVE POTENTIAL FOR A CLASS OF MODELS 
UP TO THE SECOND ORDER 

In this section, we carry out the ORSE for the system, Eq.(7), to calculate the EP up to the second order. 
The matrix elements which appear in Eq.(15) involve only Gaussian integrals except commutators of creation and 
annihilation operators and, thus, can be readily calculated as follows, 

< n\Af M [V(<f> x + *)]|m >= J2 C l n C^- n+ \n - l )!(2(2 7 r) c )-^ ±li 

v i— o 

n m n—i 

•( II f(Pi) II /(p' fc ))-ie i( ^- i+1 ^-E^- 4+I ^)^TJ^_ pi) 

j=n— i+1 k— n — i+1 2=1 

• J°° ^L e - a2 V( m - n+2i Hay/ fe 1 - h(M*) + $ + J h^J z ) (19) 

with n < m. In Eq.(19), V k (z) = d (Jjy) ■ For simplicity, in getting the above results, we have employed the 

permutation symmetry of momenta in Eq.(15) for various products of 8 functions. Note that matrix elements of (f^. 
are special cases of Eq.(19). 

Substituting the above matrix elements into Eq.(15), one can obtain the first- and the second-order corrections to 
4 O) [J;0] as 
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J X J Z 

+ J™ ^e~ a2 V(aJfJ - h(M*) + * + [ KlJz)} 



(20) 



and 



° [ ' 0J 2 J (2n)° P( P ) JJ h ** Jz 16 J (2*)° /8(p) i x 



1 



e~ ipxi h-\J z 



da 



(2n)° P[p) J X1 , 
■VW(ay/fc x 1 -I 1 (M*) + $ + J KlJ z ) 



to i!2^ y (2tt)^ nLi f(Pk)(EUi /(p*)) 



(21) 



respectively. Here, " 



represents the absolute value. 



Next, we extract the approximated EP for the system order by order. 



At the zcroth order, E^[J; cf>] = —C from Eq.(14), and so, taking — SE a^ j J -^ = h X yJ y + <f> = (p x 0) as $, one has 
= 0. Consequently, the EP at the zeroth-order of S is 



y(°)($) 



fx 



ipx—^ 



1 1 /If 2 

-f xx --Io(M 2 ) + —I 1 (M 2 ). 



(22) 



Up to the first order (hereafter, any Greek-number superscript, such as "7", "77", means "up to the order whose 
number is consistent with the Greek number"), 



El[^,S;J\ = E^[J;<f>]+5E^[J;4>] 



(23) 



and- ^"j*' M =tpi = $ yields 



7^ 



=e-*vW{ayJf£ - h(M*) + <& + / ft" 1 J,)} 



. 



(24) 



When extracting the EP up to first order, Eqs.(17), (18) and (23) imply that only the J x up to the first order, J 1 , is 
necessary. Owing to =0, it is enough to take J x — for the last term in the left hand of Eq.(24). Thus, J 1 can 
be solved from Eq.(24) as 



7-oc V 7 *" V 



(25) 



Even J 1 will not be needed to get the EP up to the first order of S, because there exists no linear, but the quadratic 
term of J x in the zeroth-order term of {Eq[J; S] — J x J x <&) as shown in Eq.(12). In fact, to obtain the EP up to the 
nth order, one need the approximated J only up to the (n — l)th order. Now one can write down the EP up to the 
first order 



■5 



V'(<&, 6) = \[f xx - /o(M 2 )] + ±M 2 h(M 2 ) - S^ 2 /" 1 

/OO J y 
- oo V 



(26) 



Obviously, this result will yield nothing but the Gaussian EP @ (1995,2002). 

Finally, we consider the second order. tp x — (p J J — — 5E o \-j^- s \ _ $ can De so l V ed for J 77 . In the present case, 
however, it is enough to use only J 7 for the EP. Substituting J 1 into Eq.(18), we obtain the EP for the system, Eq.(7), 
up to the second order as 



1 , 



v Ji (<M) = -[f xx 



I (M 2 )] + ^M 2 h(M 2 ) 



/OO J i 
^e- a2 V(aJf xx 1 ~I 1 (M 2 ) + <P) 
-oo V 71 " V 

OC J i 

^e" Q V 2 Wj- 1 - h(M 2 ) + *)] 

1 



?2 M 



16 



/i 2 - 2 



(2^ / 3 b) 



.fi fna^ 

J-oo V 71 " V 



Ii(M 2 ) + $) 



(/(E£iPfc) + ECi/(p*)) 

where, one should take S — 1 after renormalizing V 77 (<E>, 5), and /x is determined from the stationary condition 

9V 77 ($) 



. 



(27) 



(28) 



d 2 V"(<S>) 



= 



Here, V ($) is the EP after V ($, 5) is renormalized. If Eq.(28) has no real solutions, /i can be fixed by — (Qjrp 
§. Note that in (1+1) dimensions, {\[f xx - I {M 2 )} + \M 2 h{M 2 ) ~ \[i 2 f^} and [f^ - h(M 2 )} in Eq.(27) with 
S = 1 is finite and, thus, for any (l+l)-dimensional theories which make the series in Eq.(30) finite, no renormalization 
procedure is needed. 

Similarly, employing Eq.(19), one can obtain higher order corrections to the Gaussian EP from Eq.(15). To conclude 
this section, we emphasize that Eq.(27) can easily be used to give the EPs for a number of scalar field theories including 
those discussed in Refs. ||[ll} ^5|Jl5| l. 



V. APPLICATION TO \(j> 4 FIELD THEORY 

In this section, we consider the potential, 

V(cf> x ) = l -m 2 <& + \<f>i , (29) 

which was widely studied with several variational perturbation techniques |2|,|j . 

Substituting Eq.(29) into Eq.(27), and noting that a 2n+1 e~ a2 % = for n = 1, 2, • • • and ^ oo o? n e~ a ' = 
2~ n • 1 • 3 • ■ ■ (2n — 1) for n = 0, 1, 2, • • •, one can easily obtain the EP for the system, Eq.(29), up to the second order 

V 77 ($, S) = ~(f xx I (M 2 )) + hl 2 h(M 2 ) + 6(±m 2 <S> 2 + A$ 4 - \f f^) 
+S\(f xx 1 - h{M 2 ))[m 2 + 12A$ 2 + 3AC/-. 1 - Ii(M 2 ))] 



J |^7^ [m2 ^ 2 + 12A$2+ ^ 



-<5 2 12A 2 $ Z 



d u pid u P2 1 1 

(2tt)2^ f( Pl ) + /(pj) + /(p x +P2 ) /(p 1 )/(p 2 )/(p 1 + P2 ) 



s2 3 I 1 I . (30 ) 

2 j (2tt)3^ /(ELiPk)+ELi/(p*)/(ELiPfc)nLi/(p*) 
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Discarding terms with I n (M 2 )(n = 0, 1), the above result becomes identical to Eq.(2.36) in Ref. ||](1990). This can 
be verified by carrying out integrations of 7 n (0)(n = 0, 1) and I^ n \fl)(n — 2, 3, 4) in Ref. |§|(v) over one component 
of each Euclidean momentum. In (1 + 1) dimensions, taking 5—1, one finds that Eq.(30) is finite. Because the 
(I+I)-dimensional Gaussian EP in the Coleman's normal-ordering prescription is consistent M (2002) with that in 
the Stevenson's reparametrization scheme [[llj (1985), Eq.(30) with 8 — 1 and D = 1 is consistent with the Am 2 B = 
version of Sect. IV in Ref. j| (1991). As for the case of (2 + 1) dimensions, Stancu has performed a renormalization 
procedure to make V /7 (<I>) finite j| (1991). 

Using fj, fixed at the GA result for each order of the ORSE will simultaneously imply that up to each order, the 
vacuum expectation value of the field operator <j> x is identical to that in the GA. If one chooses to do so, Eq.(30) and 
also the second-order result in our former work Ref. || will yield the result in Ref. ^|(ii,iii). 

Finally, we point out that taking /i = m will lead to the conventional perturbation result on the EP. 



VI. DISCUSSION AND CONCLUSION 



In this paper, an optimized RS expansion scheme of solving the functional Schrodinger equation with an external 
source is proposed to calculate the EP beyond the GA. For the class of scalar field theories, Eq.(7), we obtain a general 
expression on the EP up to the second order which can be used easily to obtain the EP for several models. Since the 
RS expansion is a basic tool in quantum physics and the Schrodinger picture can give us some quantum-mechanical 
intuition in QFT, we believe our investigation is interesting and useful. 

Some investigations relevant to the present work can be envisioned. To our knowledge, scalar field models in 
Refs. ]l2|-|l5|] were seldom investigated beyond the GA, and so it is worth applying Eq.(27) to those concrete models 
as well as to some bosonized models in condensed matter physics. This paper discussed just ground state solution 
of Eq.(13), and actually, one can consider excited state solutions of Eq.(13) to give the EP for the excited states 
S. Our work here also implies that it is possible to introduce the optimized procedure to those schemes proposed 
in Ref. |§(i,ii,iii). Furthermore, since Ref. || has generalized the RS perturbation theory to the spinor theory, QED 
and Yang-Mills theory, it should be viable to generalize the ORSE here to those higher-spin field theories. Besides, 
since an effective action contains complete information of a field system, it will be useful to develop the ORSE here 
to calculate the effective action 0] (iii). Finally, Cornwall, Jackiw and Tomboulis developed a generalized EP for 
composite operators and calculated it with Rayleigh-Ritz procedure j^J , and it will be interesting to develop the 
ORSE here for calculating the generalized EP which will go beyond the variational result. 
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